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1 Introduction 

A theoretical understanding ol interacting lermion systems in one dimension is im- 
portant lor a number ol reasons. On the one hand, in the physics ol quasi-one- 
dimensional organic conductors [Jerome and Schulz, 1982] or ol conducting poly- 
mers [Heeger et al., 1988] interaction effects play a major role. On the other hand, 
one-dimensional models can be easier to understand than their higher-dimensional 
versions, or even exactly solvable, as is the case with the prototypical model ol corre- 
lated lermions, the Hubbard model [Lieb and Wu, 1968]. They therefore can provide 
valuable information on the role ol correlation effects in higher dimension, e.g. on the 
physics ol correlated lermions in two dimensions which is thought to be at the origin 
ol the many interesting properties ol high-temperature superconductors [Anderson, 
1987, Anderson, 1988]. 

Theoretical work on interacting lermions in one dimension has progressed along 
a number ol different lines. One approach has been the perturbative investigation ol 
the weak coupling limit. Even this is in fact not entirely straightforward, mainly be- 
cause of the infrared divergences encountered in this type of calculation which require 
a renormalization group treatment. An elementary introduction is given in chapter 2 
below, and a very complete review can be found in the literature [Solyom, 1979]. An 
alternative and more general approach is provided by the so-called "bosonization" 
method, which is based on the equivalence (valid in certain cases) between interacting 
lermions and noninteracting bosons and on the expression ol fermionic operators in 
terms ol bosons. Combined with the renormalization group approach, the bosoniza- 
tion method provides a rather straightforward description ol the peculiar properties 
ol one-dimensional interacting lermion systems ("Luttinger liquid"), and one finds 
that the low-energy physical properties are determined by only three parameters: the 
velocities ol collective charge- and spin-density oscillations (it P)CT ), and a coefficient 
K p that determines the long-distance decay ol correlation lunctions. These coef- 
ficients play a role similar to the Landau parameters ol (three-dimensional) Fermi 
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liquid theory. A number of physical properties depending on these parameters are 
discussed below, but let us mention here that in particular the coefficient K p is impor- 
tant in a much wider variety of phenomena: the temperature dependence of the NMR 
relaxation rate [Bourbonnais et al., 1984] or of X-ray scattering intensities [Pouget, 
1988], effect of impurities [Giamarchi and Schulz, 1988], or possible low-temperature 
ordered states in systems of coupled chains all depend on it. A brief discussion of 
bosonization will be given in chapter 3, but for more detailed and rigorous derivations 
and results, the reader is referred to more specialized articles [Emery, 1979, Emery, 



A rather different approach (at least until recently) is based on the famous 
"Bethe ansatz" [Bethe, 1931] which in particular has made possible an exact solution 
of both continuum fermions interacting via ^-function potentials [Gaudin, 1967, Yang, 
1967] and of the one-dimensional Hubbard model [Lieb and Wu, 1968] (and of many 
other interesting models). Without going into the mathematical details, in section 4.2 
I will try to explain the basic ingredients of the exact solution of the Hubbard model, 
and then discuss in some detail the resulting spectrum of low-lying excitations. This 
will give a rather concrete illustration of the concept of "holons" and "spinons" . 

The exact Bethe ansatz eigenfunctions are so complicated that the direct cal- 
culation of correlation functions like (3.17), (3.18) and other physical properties of 
the one-dimensional Hubbard model is hard even for very small systems 
[Ogata and Shiba, 1990] and impossible in the thermodynamic limit. In section 4.3 I 
present a method that allows in particular a determination of the coefficient K p for ar- 
bitrary correlation strength. One then obtains a rather detailed and exact description 
of the low-energy (and low temperature) properties and also of the metal-insulator 
transition occurring when the average particle number per site, n, approaches unity. 
The method generalizes rather straightforwardly to other models. 

Chapters 2 and 3 are intended to be brief and, hopefully, pedagogical introduc- 
tions to subjects where detailed reviews exist, as cited above. On the other hand, the 
material in chapter 4 is to a large extent rather recent, and in some cases previously 
unpublished. 



In this chapter we are interested in weakly interacting fermions in a one-dimensional 
metal. I consider a simple one-band model. The "kinetic energy" term in the Hamil- 
tonian (which really also contains the interaction with the static crystal potential) 
then is of the form 
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2 Weak Coupling Case 



2.1 The Model 
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where Ck )S and c\ s are the standard annihilation and creation operators for a fermion 
with momentum k and spin s. Sk is the single-particle bandstructure. In a simple 
tight -binding model one would have Sk = —2tcosk (the lattice constant is set to 
unity), but the precise form of Sk is unimportant in this chapter. The Fermi surface 
consists just of the two points ±kp. 

For weak interactions between the particles, only states in the immediate vicin- 
ity of the Fermi points are important. For these states, one then can linearize the 
electronic dispersion relation around the Fermi points, and the kinetic energy term 
takes the form 

H = v F J2{{k - M«L«M + (-k ~ k F )bl s b k , s } , (2.2) 

k,s 

Here the a (b) operators refer to states in the vicinity of -\-kp ( — kp), i.e. the a- 
particles move to the right, the 6-particles move to the left. The ^-summation is 
limited to an interval [— k 0} k ] around kp (typically, k 7r, but the precise value 
isn't important here). The Fermi velocity is given by 



VF= ^k 



(2.3) 



and the density of states per spin is N(Ep) = 1/(ttvp). 

We now introduce interactions between the fermions. For simplicity, I will only 
consider here the so-called forward (kp } s;—kp } t) — > (kp } s; — kp } t) and backward 
(kp } s; —kp } t) — > ( — kp } s; kp } t) scattering processes, with coupling constants g 2 and 
<7i, respectively. The operator decribing these processes is 

H int = — ^[gial tS bl tt a p+ 2k F + q ,th-2k F -g,s + 92ai,sK,tK+g,i a k-g,s^ ■ (2-4) 

kpqst 

Again, only processes involving particles close to the Fermi surface are important, and 
therefore the dependence of the scattering potentials </i j2 on the momentum trans- 
fer q (and possibly on the other momenta) is neglected. For Coulomb interactions 
one expects </i,</2 > 0. In principle, the long-range part of the Coulomb repulsion 
leads to a singular ^-dependence. Such singularities in g 2 can be handled rather 
straightforwardly by the techniques discussed in chapter 2, but here I shall limit 
myself to constants g\ } g2- Electron-phonon interactions can lead to effectively at- 
tractive interactions between electrons, and therefore in the following I will not make 
any restrictive assumptions about the sign of the constants. One should however 
notice that a proper treatment of the phonon dynamics and of the resulting retar- 
dation effects requires more care [Voit and Schulz, 1988]. The model thus defined 
(only forward and backward scattering, constants g\ and g 2 ) exhibits already both 
the fundamental physical phenomena of one-dimensional fermions systems, and the 
technical problems encountered in its solutions, and I will therefore concentrate on 
this model for the rest of this chapter. 
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In order to be more realistic, one can to introduce other types of interac- 
tion: on the one hand, in a half-filled band, umklapp scattering (kp, s; kp,t) — > 
( — s; — t) is possible, because the total change in momentum then is a re- 
ciprocal lattice vector. On the other hand, there is also the so-called g 4 process 
(kp, s; kp,t) — > (kp } s; kp,t). These types of interaction can lead to new phenom- 
ena in some cases, but the methods used for their treatment are similar to the ones 
discussed below for the §\ — g 2 model. 

2.2 Mean-Field Theory (beware!) 

The model introduced above exhibits various kinds of instabilities. A qualitative idea 
about this can be obtained using mean-field theory. It should however be said right 
from the outset that this is going to produce a number of qualitatively wrong results 
insofar as ground states with a broken symmetry are found, whereas exact theorems 
forbid symmetry breaking. 

To demonstrate the occurrence of instabilities in the framework of a mean- 
field theory, let us in particular consider the so-called Peierls instability, e.g. the 
instability of an interacting one-dimensional fermion gas against the formation of a 
charge density wave with wavevector 2kp- The corresponding order parameter is 

AcDW = {J2 b k,s a k+2k F ,s} ■ (2.5) 

k,s 

Note that Acdw is a complex quantity. Introducing the operator for fluctuations 
around the mean-field expactation value 

= J2 b i,s a k+2k F + q,s ~ (J2 b i,s a k+2k F + q,s) , (2.6) 
k,s k,s 

and neglecting all terms of order <5 2 , one finds 

H int ~ A C dw a l+2k F ,J 3 k,s + A* CDW ^ b\^a k+2 k F , s ■ (2.7) 

k,s k,s 

The approximate Hamiltonian then is in bilinear form, and one straightforwardly 
finds the energy eigenvalues (for k > 0) 

E k = sign(£ - k F )^/v 2 F {k-k F y + \A C Dw\ 2 , (2.8) 

and a similar expression for k < 0. One now sees that the energy of all occupied 
states is lowered if Acdw ^ 0. The total energy gain is easily calculated (for small 
Acdw) as 

AE * \Acdw\ 2 H\Acdw\/E ) , (2.9) 

and because of the logarithmic factor, it is always favorable to have |AcdvH 0- 
The most important properties of the state with |AcdvH are: 
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1. The average charge density is modulated, i.e. one has a charge-density wave 
(CDW): 

(p(x)) = po + Sp cos(2k F x + <j>) , (2.10) 

where <f> is the phase of Acdw, and p is the average charge density. The CDW 
thus break, s the continuous translational symmetry of the §\ — g 2 model, vari- 
ations of <f> representing continuous translations of the CDW. The breaking of 
the continuous symmetry implies that there are low-energy collective excita- 
tions (Goldstone modes) with a linear energy-momentum relation u(k) = c\k\, 
with some velocity c. These excitations represent long-wavelength fluctuations 
of the wavelength of the CDW. 

2. There is a gap in the spectrum of single particle excitations. Consequently, 
there is no Pauli susceptibility, and no single-particle conductivity. However, 
the CDW can move as a whole, leading to the charge-density wave conductivity 
first proposed by Frohlich. 

One can investigate other types of instability along similar lines. The relevant 
operators which can take on non-zero expectation values are: 

1. spin density wave (SDW): 

S DW a = ^2 bl )S cr" tt a k+2 k F ,t , (2-H) 

k,s,t 

2. singlet superconductivity (SS): 

OSS = J2 sb -k -sdk,s , (2.12) 
k,s 

3. triplet superconductivity (TS): 

T s a = ^2sb_ k _ s (T^ t a k ^ t , (2.13) 

k,s,t 

where o a are the Pauli spin matrices (a = x } y } z). These operators are Fourier 
transforms of the real space operators 

OcDw(x) = ^2%l) ] _^(x)^+,s{x) , 

S 

S DW a (x) = J^^-A^^t^+Ax) > 

s,t 

Oss(x) = J2 st l ; --s( x ) t l ; +A x ) i 

s 

OtsJx) = Y.^-M. x XtM x ) , (2-14) 
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Figure 1: Ground state phase diagram of the g\ — g 2 model obtained in the mean-field 
approximation 

where ip TtS (x) is the annihilation operator for a right-going (r = +) or left-going 
(r = — ) particle with spin s. 

Seeking the energetically most stable state in the g\ — g 2 ~plane, one finds the 
"phase diagram" shown in fig. f . All the states in fig. f break continuous symmetries: 
in the CDW and SDW states, translational symmetry is broken, in the SS and TS 
states gauge symmetry is broken. In addition, in the "triplet" phase SDW and TS, 
spin rotation invariance is also broken. Consequently, all theses states have at least 
one Goldstone mode in their excitation spectrum. 

Compared to a three-dimensional interacting fermion system, fig. f exhibits a 
remarkable feature: in no part of the parameter space (except in the trivial singular 
point gi = g 2 = 0) is a "normal" (Fermi liquid like) state stable. While the pairing 
instabilities SS and TS are indeed the analogue of the well-known Cooper insta- 
bility occuring in any dimension for weakly attractive forces between fermions, the 
density-wave type instabilities (CDW, SDW) are due to the particular form of the 
one-dimensional Fermi "surface" , and have no direct analogues in higher dimensions 
(apart from special cases with nesting Fermi surfaces). In particular, weakly repul- 
sive fermions in three dimensions are Fermi liquids without any signs of symmetry 
breaking. 

One can extend the mean-held calculation in order to check its internal con- 
sistency. A particular useful check is a calculation of the fluctuations of the order 
parameter around its mean value. Without going into details, one finds 



Here u q = c\q\ is the energy of the Goldstone modes, f3 = 1/T, 7 = CDW, SDW, 
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SS, or TS, and £ is the BCS coherence length for the pairing instabilities or an 
analogous length for the density wave instabilities. Expression (2.15) is divergent 
in one dimension both at zero and at finite temperature (and in two dimensions for 
T / 0), indicating that in these cases fluctuations of the order parameter about 
its supposed mean value are infinitely large. This clearly indicates that something 
went wrong in the mean-held theory, and that the assumption of a symmetry-broken 
ground state doesn't make sense. A more consistent treatment will be given in the 
following sections. 



2.3 Renormalization Group: Coupling Constants 

In order to understand the problems with the mean-held approach discussed at the 
end of the last section, it is helpful to take a step back and to look at the usual micro- 
scopic justification of the instability criteria used. The typical case is the BCS theory 
of superconductivity (in three dimensions). Here the mean-held critical temperature 
corresponds to a divergence of the pair susceptibility 

A(T) = J drdT{T T 4>l(x, t)^_ s (x, t)V> s (0, 0)V>- s (0, 0)) (2.16) 

calculated in the random phase approximation. One then has 

A(T) = Ao(T)/(l + AA (T)) , (2.17) 

where A (T) is the noninteracting pair susceptibility. The standard justification for 
the approximation (2.17) lies in the fact that A (T) is logarithmically divergent at 
low temperatures, and that the RPA takes into account the most divergent diagrams 
at each order in perturbation theory. 

In one dimension, things are more complicated: in addition to the logarithmic 
divergence in A (T), the density-density correlation function at wavevector 2kp 

n (T) = j J2 /(£fc) ~ /(£fc " 2%) (2.18) 

L k £fc — £k-2k F 

is also logarithmically divergent. In a diagrammatic language, A and il are given 
by the particle-particle and particle-hole diagram in fig. 2. Consequently, at each 




Figure 2: The elementary logarithmically divergent particle-particle and particle-hole 
diagrams 

order in perturbation theory, a whole set of mixed particle-particle and particle- 
hole diagrams has to be summed, leading to the so-called parquet equations. In the 
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present context, this has been done by Bychkov et al. in 1966 [Bychkov et al., 1966]. 
However, instead of discussing their method, I will here discuss a renormalization 
group approach which is based on Anderson's "poor man's scaling" of the Kondo 
problem (where similar mixed divergences occur) [Anderson, 1970, Anderson et al., 
1970]. This approach has been adapted to the one-dimensional fermion problem by 
Solyom [Solyom, 1979], and I here closely follow his presentation. 

The basic idea is by now well known: the divergent terms are of the form 
ln(E /T) (or ln(E /io) } ln(E / vp\q\) at T = 0), where E = vpk . Consequently, if it 
is possible, via an elimination of high-energy degrees of freedom, to find a model with 
a smaller cutoff energy E, but the same physical properties at energies much lower 
than E, the problem of divergent diagrams is less serious. Moreover the problem may 
even be completely cured if it is possible to apply the transformation repeatedly until 
E ~ T (provided one still remains within the range of applicability of perturbation 
theory). 

Physical properties at low energies are determined by the two-particle T-matrix 
T(lo) = H mt + — T (2.19) 

U> + 16 — ti 

which is related to the transition probability between an initial state \i) and a final 
state \f) via 

W i ^ f = 2ic\(i\T(E i )\f)\ 2 p(f) , (2.20) 

where E{ is the energy of the initial (and final) state, and p(f) is the density of final 
states. Consider initial and final states 

io = <A,/*i > i\f) = b LAJ°) > (2-2i) 

where |0) is the filled Fermi sea, and £4, & 4 kp } k 2} k 3 — kp. Second order pertur- 
bation theory in Hi nt gives 




a,8<Jf3r/ 



(2.22) 



where in the summation over intermediate states of energy E only the contributions 
divergent for u — > have been retained. Note that the only nonzero contributions 
come from the ^-interactions, all singular contributions from g 2 cancel each other. 

The perturbative result (2.22) can now be used to generate renormalisation 
group equations. We want to reduce the cutoff E to a lower value Eq } and simul- 
taneously adjust the coupling constants with the requirement that the physics at 
low energies u <^ E' remains unchanged, i.e. that 

T^7(^HI/) = . (2.23) 
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Using the explicit result (2.22) and changing the name of the upper cutoff E' — >■ E), 
one finds the renormalization group equations 

E^9i{E) = —gl(E) , 

(IE TTVp 

E-^92(E) = -^gl(E) . (2.24) 
(IE Zttvf 

The initial conditions for these equations are that when the "running" cutoff E equals 
the initial cutoff E 0} the coupling constants should take their bare values: gi(E ) = g^. 

One can look at this derivation from a slightly different point of view: per- 
turbatively, the gi-like term in (i\T(uj)\f) (i.e. the term proportional to 8 an 8pj) 
is 

T 1 {u)=g 1 + ^-\n{E lu) , (2.25) 



which can be rewritten, in terms of a reduced cutoff energy 



g 1 + jL HEo/E ' o] 

TTVp 



+ jLi n ( E 'Ju>) . (2.26) 

TTVp 



This last expression now can be interpreted as a perturbative result for a system 
where the initial coupling constant and cutoff (</i, E ) are replaced by new parameters 
{g'^K), with 

g[ =9i + J ^HE /E' ) . (2.27) 

Using an E' infinitesimally smaller than E 0} and iterating the transformation (</i, E ) 
— > ((/^-Eq), one recovers the renormalization group equations (2.24). 

The renormalization group equations (2.24) show that there is no renormaliza- 
tion of the coupling constants if g\ = 0. The general solution for g\ is 

= i + aw) ' (2 ' 28) 

and g2(E) is obtained straightforwardly noting that the combination g\ — 2g 2 is un- 
renormalized. This result now shows that there are two qualitatively different regimes: 
if initially g\ > 0, one has g\(E) — > when the cutoff energy decreases, i.e. one 
renormalizes to a weak-coupling regime. One then expects that the perturbative 
approach used here should give reliable results for the low-energy properties of the 
model. On the other hand, if initially g\ < 0, gi(E) diverges at some finite cutoff 
energy. Clearly, close to the divergence, the assumption of weak coupling underlying 
the present expansion breaks down. The physical properties in this regime can how- 
ever be understood using the bosonization method introduced in the next chapter. 
The cases g\ > and g\ < correspond to the ferromagnetic and antiferromagnetic 
Kondo problems, respectively. 
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2.4 Renormalization Group: Correlation functions 

Let us now use the above results for §\ > 0, where everything is well-defined, to calcu- 
late correlation functions. For dehniteness, let us consider the correlation function of 
charge-density fluctuations with wavevector around 2kp, as represented diagrammat- 
ically in fig. 3. In RPA, this function diverges at the mean-held critical temperature 








Figure 3: Perturbativ expansion of the CDW correlation function T>cdw (the hatched 
bubble) to first order in the §{. 



of the Peierls instability. To zeroth order one finds, at T = 

<n° < \-rP( \ 1 i l" 2 + "h 2 ] 1 ' 2 



TTVp 

The hrst order result is 

Vcdw = T>°[1 + {g C Dw/2)V° + ...] 



E 



(2.29) 



(2.30) 



with gcDW = 2(/i — §2- For the other types of instabilities (SDW, SS, TS), analogous 
results can be found, with coupling constants gsDW = —92^9ss = 9\ + 92^9ts = 
-9i +92- 

We now use these perturbative results to obtain renormalization group equa- 
tions for the correlation functions Z\. The general structure of the perturbation 
theory is, at energy cutoff E 



1 ln(uj/E)[l - \n(u/E) + ...] 



P -- 

TTVp 2tTVf 

and therefore varying the cutoff, one has 
irvpdVi = [-1 + —ln(u/E)]^ r . 

TTVp E 



12.31] 



(2.32) 



We now have to take into account two effects: (i) if the running cutoff E is smaller 
than the initial cutoff E 0} a certain number opf states have already been eliminated, 
and therefoe the coupling constant §i appearing in (2.32) actually is the renormalized 
coupling constant, (ii) for the same reason, the factor f /(tvvf) \n(u/E) appearing in 



fO 



(2.32), which in a purely perturbative calculation represents Z>°, should be replaced 
by the full correlation function at cutoff E. The renormalization group equation for 
Di then becomes 

dV- 

™ f !)^e = ~ 1 + 9 ^ V ^ ■ ^ 

An initial condition is provided by the observation that when the cutoff is of the 
order of the external frequency u, the correlation function are finite and small, i.e. 
T>i(io; E io) (the precise value is of minor importance here). One then can 
integrate eq. (2.33), starting with the initial condition until the cutoff reaches E . 
The leading singular behaviour at low energies and momenta is found to be 

Vi(q,u) « [to 2 + v 2 F q 2 ]- a < /2 ^'([u 2 + v 2 F q 2 ]/E 2 ) . (2.34) 

Here the exponents V{ are 

CtCDW = CtSDW = (#2 - gi/2)/(TTV F ) , 

ass = a TS = -(g 2 - gi/2)/(Trv F ) , 
VSDW = vts = 1/2 , 
vcdw = vss = —3/2 . 

The important result here is the existence of power laws in the low-energy, low 
momentum regime, with the value of the exponents depending on the interaction 
constants, i.e. the exponents are non-universal. If a 8 - > 0, the corresponding cor- 
relation function is enhanced over the logarithmic behaviour of the noninteracting 
case. Moreover, because of the logarithmic correction factors in (2.34), the SDW and 
TS correlations are always enhanced over the CDW and SS. The "phase diagram", 
indicating the strongest fluctuations in different parts of the §\ — g 2 plane, is shown 
in fig. 4. 

The real space correlation functions are found Fourier transforming (2.34) as 

Vi(r,T) w (r 2 + 4r 2 )-( 1 - ai/2 > In"' (r 2 + v 2 f t 2 ) , (2.36) 

i.e. as expected there is no long range order, even at zero temperature. However, pro- 
vided that a.i > 0, the spatial decay of correlations is slower than in the noninteracting 
case. Similary, at finite temperatures, one finds Z\-(0,0) ~ (T/E )~ a * \n v ' (E /T) } i.e. 
there are no singularities, and therefore no phase transitions at finite temperature, 
but the fluctuations with the strongest exponent a 8 - dominate at low temperature. 
Contrary to the correlation functions Z\, thermodynamic properties like the spin sus- 
ceptibility, the specific heat, or the compressibility have only finite corrections at low 
temperatures. 

The calculations presented here are based on lowest nontrivial order pertur- 
bative calculations, and consequently, the exponents a 8 - are correct to order §{. It 
is possible to go to higher order in the calculation. One then faces an interesting 
new phenomenon: self-energy corrections to the one-particle Green's function be- 
come singular at low energies and momenta (they are regular at first order). This 
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Figure 4: Zero temperature "phase diagram" showing the most divergent type of flutu- 
ations in the plane. Phases in parentheses are logarithmically weaker divergent. 

signals the breakdown of the Landau quasiparticle picture in one dimension: in fact, 
as discussed in the next chapter, the elementary excitations at low energy are col- 
lective charge and spin fluctuations, rather than individual quasiparticles. As far as 
the renormalization group calculation is concerned, this means that it is no longer 
straightforward to define transition amplitudes between initial and final quasiparti- 
cle states as in (2.21, 2.22), because there are singular corrections to these states 
themselves. This situation can be handled by a field-theoretical formulation of the 
renormalization group, as described in detail by Solyom [Solyom, 1979] . This type of 
approach has been pursued up to third order in the coupling constants. One should 
however notice that in higher order, some of the integrations in the diagrammatic 
calculation are nonsingular, and therefore necessarily depend on details of the elec- 
tronic bandstructure and on the cutoff procedure used. It is thus not very promissing 
to push this approach to higher and higher order if one wants to understand strongly 
interacting fermions. On the other hand, the bosonization method allows one to ob- 
tain a number of general exact statements, which in some cases (e.g. the Hubbard 
model, see below) can be used to extract all the parameters determining the low 
energy physics, even in strongly interacting cases. 
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3 Bosonization, spin— charge separation, 
Luttinger liquid 

3.1 Bosonization Formalism 

One of the important results in the theory of one-dimensional interacting fermions 
is that fermion field operators can be expressed in terms of boson operators 
[Luther and Peschel, 1974, Mattis, 1974]. This equivalence then can be used to ex- 
press the fermion Hamiltonian in a particularly simple form, in terms of boson fields 
only (see eq.(3.3) below). Consider for the moment spinless electrons, and define 
a boson field <f>(x) by d(f>/dx = Trp(x) } where p is the deviation from the average 
density. Then, creating a particle at site x means introducing a kink of height tt 
in </>, i.e. at points on the left of x cf> has to be shifted by tt. Displacement opera- 
tors are exponentials of momentum operators, and therefore a first guess would be 
ip^x) exp(z7r j^oo H(x')dx'), where II is the momentum density conjugate conjugate 
to <f>: [<f)(x),H(y)] = iS(x — y). However, this operator commutes with itself, instead 
of satisfying canonical anticommutation relations. Anticommutation is achieved by 
multiplying with an operator, acting at site x, that changes sign each time a particle 
passes through x. Such an operator is exp(H<f)(x)). The final result then is 



t -1 r . f x 

ip±(x) = lim — -j= exp ±ikpx^fi(f>(x) + iir / Tl(x')dx' 



(3.1; 



where the upper and lower sign refer to electrons near kp and —kp respectively, and 
a is a short-distance cutoff. A detailed derivation of this important relation as an 
operator identity is given in the literature [Haldane, 1981a, Heidenreich et al., 1980]. 
One should also notice that (3.1) is strictly speaking valid for models with linear 
energy-momentum relation. If there is some curvature in the dispersion relation (as 
is necessarily the case in lattice models), all odd powers of exp(z</>) appear [Haldane, 
1981b], i.e. 



v4 



2^ a m exp 

m>0 



±z(2m + l)(k F x - <f>(x)) + nr J\{x')dx' 



(3.2) 



where a m are model dependent constants. 

For electrons with spin, one simply introduces one boson field for each spin 
orientation. Introducing charge and spin bosons via </> Pj(7 = (</>j ± <f>^)/\/2, one then 
finds that the low-energy, large-distance behaviour of a one-dimensional fermion sys- 
tem with spin-independent interactions is described by the Hamiltonian [Solyom, 
1979, Emery, 1979] 

H = H P + H a + — ^- / dx cos(V8^) . (3.3) 

Z7TO! J 
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Here a is a short-distance cutoff, g\ is the backward scattering amplitude, and for 
v = p,a 

'ttu,,K,,~ u,. 



H,, = dx 



2 v 2ttK„ 
The full expressions for the phase fields are* 



(3.4) 



M*) = -^^-^(pj + ./.Wl-JV^, (3.5) 
L '~j n p L 

n„(x) = We- alpH2 - tp nMp)-v-(p)] + JJL . (3.6) 

L p*0 

Here p r (p) (cr r (p)) are the Fourier components of the charge (spin) density operator 
for right-(r = +) and left-(r = — ) going fermions. Introducing the total number 
operators (measured with respect to the ground state) N rs for right- and left-going 
particles (r = ±) of spin s, the (charge and spin) number and current operators are 
N v = [{N + ^N^)±{N +i + N_ i )]I^J v = [{N + ^-N^)±{N +i -N_ i )]l^ where 
the upper and lower sign refer to charge and spin, respectively. 

The operators cf> v and in (3.3) obey Bose-like commutation relations: 

[4> v {x),U ll {y)] = i8 vll 8{x-y) , 

and consequently, at least for g\ = 0, (3.3) describes independent long-wavelength os- 
cillations of the charge and spin density, with linear dispersion relation u v (k) = u v \k\, 
and the system is conducting. This model (no backscattering), usually called the 
Tomonaga-Luttinger model, is the one to which the bosonization method was origi- 
nally applied [Luttinger, 1963, Mattis and Lieb, 1965]. The only nontrivial interac- 
tion effects in (3.3) come from the cosine term. However, for repulsive interactions 
(gi > 0), this term is renormalized to zero in the long-wavelength limit, and at 
the fixed point one has K* = 1. The three remaining parameters in (3.3) then 
completely determine the long-distance and low-energy properties of the system. 
It should be emphasized that (3.3) can be derived exactly for fermions with linear 
energy-momentum relation. For more general (e.g. lattice) models, (3.3) is still the 
correct effective Hamiltonian for low-energy excitations. 

The Hamiltonian (3.3) also provides an explanation for the physics of the case of 
negative where the renormalization group scales to strong coupling (cf. sec. 2. 3). 
In fact, if \gi\ is large in (3.3), it is quite clear that the elementary excitations of 
cf><j will be small oscillations around one of the minima of the cos term, or possibly 
soliton-like objects where cf> a goes from one of the minima to the other. Both types of 
excitations have a gap, i.e. for §\ < one has a gap in the spin excitation spectrum, 
whereas the charge excitations remain massless. This can actually investigated in 
more detail in an exactly solvable case [Luther and Emery, 1974]. 



"The N v - and J„-terms are discussed by Haldane [Haldane, 1981a] 
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In a half-filled band, umklapp scattering is possible, transferring two particles 
from — kp to kp, or vice versa. These processes lead to an extra term 

/0 2g \ 9 / dx cos(y/S(f) p ) (3.7) 

in the Hamiltonian. Similarly to the g\ term, this can lead to a gap, this time in 
the charge excitation spectrum [Emery et al., 1976, Solyom, 1979]. The ground state 
then is insulating. This happens in particular for the Hubbard model with repulsive 
interaction at exactly half-filling. 



3.2 Spin— Charge Separation 

One of the more spectacular consequences of the Hamiltonian (3.3) is the complete 
separation of the dynamics of the spin and charge degrees of freedom. For example, in 
general one has u a ^ u p , i.e. the charge and spin oscillations propagate with different 
velocities. Only in a noninteracting system one has u a = u p = vp. To make the 
meaning of this fact more transparent, let us create an extra particle in the ground 
state, at t = and spatial coordinate x . The charge and spin densities then are 
easily found, using p(x) = — (\/2/Tr)d(f> p /dx (note that p(x) is the deviation of the 
density from its average value) and a z (x) = — (\/2 / %)d(f) a / dx : 

(0\tjj + (x o )p(x)'ijj\_(x o )\0) = S(x - x ) , 
(0|^(zoK(z)V>i(z o )|0) = 6{x-x ) . (3.8) 

Now, consider the time developement of the charge and spin distributions. The 
time-dependence of the charge and spin density operators is easily obtained from 
(3.3) (using the fixed point value g\ = 0), and one obtains 

(0\tlj + (x o )p(x,t)tlj\_(x o )\0) = S(x-x -u p t) , 
(0\tl) + (x o )(j z (x,t)tlj\_(x o )\0) = S(x-x -u a t) . (3.9) 

Because in general u a ^ u p , after some time charge and spin will be localized at 
completely different points in space, i.e. charge and spin have separated completely. 
A interpretation of this surprising phenomenon in terms of the Hubbard model will be 
given in sec. (4). Here a linear energy-momentum relation has been assumed for the 
electrons, and consequently the shape of the charge and spin distributions is time- 
independent. If the energy-momentum relation has some curvature (as is necessarily 
the case in lattice systems) the distributions will widen with time. However this 
widening is proportional to -v/t, and therefore much smaller than the distance between 
charge and spin. Thus, the qualitative picture of spin-charge separation is unchanged. 



3.3 Luttinger Liquid 

The simple form of the Hamiltonian (3.3) at the fixed point = makes the calcula- 
tion of physical properties rather straightforward. First note that acoustic phonons in 
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one dimension have a linear specific heat. Consequently, the low-temperature specific 
heat of interacting fermions is C(T) = ~{T } with 

7/70 = t^( v f/u p + v F /u a ) . (3.10) 

Here 70 is the specific heat coefficient of noninteracting electrons of Fermi velocity 
v F . 

The spin susceptibility and the compressibility are equally easy to obtain. Note 
that in (3.3) the coefficient u a / ' K a determines the energy necessary to create a nonzero 
spin polarization, and u p /K p fixes the energy needed to change the particle density. 
Given the fixed point value K* = 1, one finds 

x/Xo = v F /u rj , k/kq = v F K p /u p , (3-11) 

where Xo an d K o are the susceptibility and compressibility of the noninteracting case. 
From eqs.(3.10) and (3.11) the Wilson ratio is 

R w = ^ = p — . 3.12 

7 XO Up + Ua 

The quantity Tl p (x) is proportional to the current density. Obviously, the Hamil- 
tonian commutes with the total current, and therefore the frequency dependent con- 
ductivity is a delta function at u = 0. Using the Kubo formula, one straightforwardly 
finds 

a(u) = 2K p u p 6(u) , (3.13) 

i.e. the product K p u p determines the weight of the dc peak in the conductivity. 

The above properties, linear specific heat, finite spin susceptibility, and dc con- 
ductivity are those of an ordinary Fermi liquid, the coefficients and K p de- 
termining renormalizations with respect to noninteracting quantities. We will now 
consider quantities which show that a one-dimensional interacting fermion system 
is not a Fermi liquid. Consider the single-particle Green's function which can be 
calculated using the representation (3.1) of fermion operators. One then finds for the 
momentum distribution function in the vicinity of k F : 

n k n kp - /3sign(& - k F )\k - k F \ a , (3-14) 

and for the single-particle density of states: N(lo) |^|", with a = (K p -\-l/ K p — 2)/4, 
and f3 is a model-dependent constant. Note that for any K p ^ 1, i.e. for any 
nonvanishing interaction, the momentum distribution function and the density of 
states have power-law singularities at the Fermi level, with a vanishing single particle 
density of states at E F . This behaviour is obviously quite different from a standard 
Fermi liquid which would have a finite density of states and a step-like singularity 
in rik- The absence of a step at k F in the momentum distribution function implies 
the absence of a quasiparticle pole in the one-particle Green's function (this can of 
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course be verified directly). The unusual type of behaviour found here has been called 
Luttinger liquid by Haldane [Haldane, 1981a]. 

The coefficient K p also determines the long-distance decay of all other corre- 
lation functions of the system: Using the representation (3.1) the operators (2.14) 
become 



OcDwix) = lim- e-*^<^cos[^(x)] , (3.15) 

OsdwM = lim^-e-^^ cos[V2e„(x)] , (3.16) 

where 9 v (x) = tt J x H l/ (x')dx l . Similar relations are also found for the other operators 
in (2.14). It is important to note here that all these operators decompose into a prod- 
uct of one factor depending on the charge variable only by another factor depending 
only on the spin field. Using the Hamiltonian (3.3) at the fixed point gl = one 
finds for example for the charge and spin correlation functions^ 

(n(x)n(O)) = Kp/iirx) 2 + A 1 cos(2k F x)x- 1 - Kp ln- 3/2 (x) 

+ A 2 cos(4:k F x)x- 4Kp + ... , (3.17) 
{S(x)-S(0)) = 1/(ttx) 2 + B 1 cos(2k F x)x- 1 - Kp ln 1/2 (x) + ... , (3.18) 

with model dependent constants Ai 7 B{. The ellipses in (3.17) and (3.18) indicate 
higher harmonics of cos(2k F x) which are present but decay faster than the terms 
shown here. Similarly, correlation functions for singlet (SS) and triplet (TS) super- 
conducting pairing are 

(Ol s (x)Oss(0)) = Cx- 1 - 1 ^ln- 3 / 2 (x) + ... , 
(OU^OtsM) = Dx- 1 - 1 / Kp ln 1 / 2 (x) + ... . (3.19) 

The logarithmic corrections in thses functions have been found rather recently 
[Giamarchi and Schulz, 1989]. The corresponding susceptibilities (i.e. the Fourier 
transforms of the above correlation functions) behave at low temperatures as 

Xcdw(T) « T^-'IMT)^ 2 , X sdw(T) « T K ^\ ln(T)! 1 / 2 , (3.20) 

Xss(T) « T^-'lHT)^ 2 , X ts(T) « T^-'IMT)^ 2 , (3.21) 

i.e. for K p < 1 (spin or charge) density fluctuations at 2kp are enhanced and diverge 
at low temperatures, whereas for K p > 1 pairing fluctuations dominate. The remark- 
able fact in all the above results uis that there is only one coefficient, K p , which 
determines all the asymptotic power laws. 

If in a half-filled band umklapp scattering leads to an insulating state with a 
gap in the charge excitation spectrum, the above results are valid with K p = (then 

indicates exponential decay). 



^The time- and temperature dependence is also easily obtained, see [Emery, 1979]. 
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4 The Hubbard model in one dimension 
4.1 The Hamiltonian and its symmetries 

The Hubbard model is the protypical model used for the description of correlated 
fermions in a large variety of circumstances, ranging from high-T c superconductors 
to heavy fermion compounds and organic conductors. In spite of its apparent sim- 
plicity, there is still no general solution, or even a consensus on its fundamental 
properties. Notable exceptions are the cases of one and infinite dimensions [Voll- 
hardt, 1992, Georges and Kotliar, 1992]. In particular, in one dimension an exact 
solution is available. This solution gives exact energies of the ground state and all 
the excited states in terms of the solution of a system of coupled nonlinear equations. 
On the other hand, the corresponding wavefunctions have a form so complicated that 
the explicit calculation of matrix elements, correlation functions and other physical 
quantities has remained impossible so far. In the following sections I shall describe 
the features of the one-dimensional Hubbard model which make it exactly solvable, 
and subsequently describe in detail the energy spectrum. Finally, I will show how 
the knowledge of the energy spectrum can be combined with the results of the pre- 
ceding two chapters to obtain a rather detailed picture of the low-energy properties, 
in particular of correlation functions. 

The Hamiltonian in one dimension has the well-known form 



where c 8jS is the fermion annihilation operator on site i with spin s, n 8jS is the cor- 
responding number operator, and the sum is over the L sites of a one-dimensional 
chain with periodic boundary conditions. 

The model has two global SU(2) symmetries [Pernici, 1990, Zhang, 1990, Schulz, 
1990a]: the first is the well-known spin rotation invariance, with generators 



The second type of symmetry is particular to the Hubbard model and relates sectors 
of different particle numbers. Its generators are 



The total symmetry thus is SU{2) X SU{2) ~ 5'(9(4). One should notice that more 
complicated interactions, e.g. involving further neighbours, will conserve the spin 
rotation invariance but in general not the "charge" SU(2) invariance (4.3). Rather, 
this second symmetry will become the standard global U(l) invariance associated 
with particle number conservation. It is nevertheless possible to construct particular 
types of further-neightbour interactions which do conserve the full SU (2) X SU (2) 
invariance. 








(4.3) 
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4.2 The exact solution: a brief introduction 

I here closely follow a presentation by N. Andrei [Andrei, 1988]. To obtain the 
wavefunctions of the Hubbard model, one writes a general A^-particle state as 

N 

\ F ) = E • • • E E • • • E f^-,A^ ' x ■ II c L» \ vac ) ■ ( 4 - 4 ) 

Xl X N S! S N i = \ 

The simplest case is the two-particle problem. Then in the two parts of configuration 
space x 2 > Xi (region I) and X\ > x 2 (region II) the wavefunction is just a product 
of two plane waves, and the only nontrivial effects occur for X\ = x 2 . The full 
wavefunction F then is 

Fs uS2 = e^ +k ^[6(x 2 - Xl )H uS2 + 6( Xl - x 2 )Hlj > (4-5) 

where the factors £ depend only on the spin quantum numbers of the two particles. In 
order to satisfy the Schrodinger equation, these coefficients have to obey the equation 

^ll^ = E ^s lS2 ,t t t 2 C I t 2 • (4-6) 
tl,t 2 

The S-matrix has the form 

_ t(sm k x - sin k 2 ) + ill P a p :1 s 
i(sin ki — sm k 2 ) + lU 

where P a f5 ri s is the operator permuting two spins. This operator acts on the part of 
the two-particle Hilbert space. 

For N particles, there are AH different regions in configuration space, each 
correspond to a permutation Q of the sequence (1, 2, 3, ... , N), so that xq\ < xq 2 < 
. . . < xqn. The wavefunction takes the form 

F Su ..., SN (x 1 ,...,x N ) = e i Zi k >*> OqMU,..,sAQ) > (4-8) 

Qes N 

where Sn is the group of permutations of A^ objects, and = 1 if xq\ < xq 2 < 

. . . < xqn and zero otherwise.. The Schrodinger equation imposes the condition 

Z(Q') = s t3 C(Q) (4-9) 

when Q and Q' are two permutations which only differ by the exchange 

However, for more than two particles, there is more than one way to go from 
one part of configuration space to another. Consider three particles. Denoting by 
(ijk) the part of configuration space with X{ < Xj < Xk, there are obviously the 
possibilities 

(123) -> (213) -»■ (231) -»■ (321) 

and 

(123) -> (132) -»■ (312) -»■ (321) 
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In order for the two paths to lead to the same amplitude in (321), one has to have 



S™S"S 12 = S L2 S Ld S™ . (4.10) 

These are the famous Yang-Baxter equations [Yang, 1967] which have to be satisfied 
for a system to be solvable by Bethe Ansatz. One can verify that these equations are 
indeed satisfied by the S-matrix of the Hubbard model, eq. (4.7). 

One then imposes periodic boundary conditions. This means that the relation 

F Sl ,...,S N ( X l 7 • • • 7 X j 0, . . .) F Sl ,...,S N ( X l 7 • • • 7 X j L, . . .) (4.11) 

should be valid for every j. When Xj varies from to L two things happen: (i) the 
wavefunction is multiplied by a factor e lklL , and (ii) particle j is permuted with all 
the other particles. To satisfy periodic boundary conditions, one therefore needs to 
satisfy the relation 

Z 3 U,.,s N (D = ^ L U...,s N (I) , (4-12) 

where Zj, acting on the spin Hilbert space of N particles, is an operator representing 
the phase shifts as particle j crosses all the other particles: 

Zj = S 3 ' 3 ' 1 . . . S^S''" . . . S 3 ' 3+1 (4.13) 

Thus, the spin wavefunctions £ are the eigenfunctions of the operator Zj, and the 
allowed values of the momenta k are related to the eigenvalues of Zj. The deter- 
mination of the eigenfunctions of Zj is related to the so-called 6-vertex model of 
statistical mechanics. A detailed description of the method used can be found in 
more specialized articles [Baxter, 1982, Thacker, 1982]. 

Using those results, the allowed values of kj are obtained from the solution of 
the coupled set of nonlinear equations 

e ik 3 L = || e ( -v---^ | ( 414 ) 

(4.15) 

Here N is the total number of electrons, M is the number of down-spin electrons 
(M < N/2), and e(x) = (x + i)/(x — i). The X a are parameters characterizing the 
spin dynamics. We note that in general, both the &j's and the A's are allowed to be 
complex. The energy and momentum of a state are 

N N 

E = -2tJ2cosk 3 , P = J2k 3 . (4.16) 

j=i j=i 
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4.2.1 Solutions of the Bethe ansatz equations 



The determination of all the solutions of eqs. (4.14, 4.15) is not easy. It has recently 
been shown (under certain assumptions) that these equations do indeed give all the 
"lowest weight" (with respect to SU{2) X SU{2)) eigenstates of the Hubbard model, 
i.e. all states satisfying r)\ip) = Cl^} = 0- The complete set of eigenstates then is 
obtained acting repeatedly with or C on |V>) [Essler et al., 1991]. H ere I will 
limit myself to the ground state and to the low-lying elementary excitations. These 
questions have been investigated in some detail [Shiba, 1972, Coll, 1974, Woynarovich, 
1983], however, the finite chain data presented below seem to be quite useful in 
understanding the nature of the excitations, and I therefore discuss them in detail. 

If both the &'s and the A's are all real, only the phases in (4.14, 4.15) have to 
be determined. Taking the logarithm of these equations, one finds 

M 

L]z 3 = 2kI 3 + 2 J2 arctan[4(A„ - sin kj)/U] (4.17) 

a = l 

N M 

2^arctan[4(A„ - sin %)/[/] = 2ttJ„ + 2 ^ arctan[2(A„ - Xp)lU] . (4.18) 

3=1 f3=l 

The quantum numbers {Ij} are all distinct from each other and are integers if M is 
even and half-odd integers (HOI, i.e. of the form 1/2, 3/2, . . .) if M is odd, and 
are only defined modulo L. Similarly, the set { J a } are all distinct and are integers if 
N — M is odd and HOI if N — M is even. Moreover, there is the restriction 

\J a \ < (N-M + l)/2 . (4.19) 

Summing (4.17) over j and (4.18) over a, the total momentum is found as 

p = Y + • (4-20) 



Ground state. The ground state is nondegenerate only if N is of the form 4z/ + 2 (y 
integer): obviously, if N is odd, the ground state has (at least) spin 1/2. Further, if 
N is an integer multiple of 4 the noninteracting ground state has a sixfold degeneracy, 
and in the interacting case the ground state turns out to be a spin triplet. In the 
following I shall restrict myself to the case of N = N = 4z/ + 2, i.e. the ground 
state is nondegenerate (in the following, N will denote the particle number in the 
ground state). The ground state then is a singlet, with M = iV /2, i.e. M is odd. 
The allowed values of the J's range from —(N /2 — l)/2 to (N /2 — l)/2. There are 
exactly N /2 such integers, i.e. all the J's are fixed. The J's are consecutive between 
-(N - l)/2 and (N - l)/2, i.e. 

ih) = {-(iV -l)/2,...,(iVo-l)/2} , (4.21) 
{J a } = {-(iV /2-l)/2,...,(iVo/2-l)/2} . (4.22) 
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In the thermodynamic limit L — > oo the distance between consecutive &'s or A's 
decrease like 1/L, and one can then find linear integral equations for the density 
of &'s and A's on the real axis. Numerical results for the ground state energy as a 
function of particle density and U have been given by Shiba [Shiba, 1972]. 

"4kF" singlet states. Excited states are obtained by varying the quantum num- 
bers. The first possibility, giving rise to excited singlet states, is obtained by removing 
one of the 7's from the ground state sequence (4.21) and adding a "new" 7: 

ih) = {-(iV -l)/2,...,-(iVo-l)/2 + zo-l, 

_(iV -l)/2 + zo + l,...,(iVo-l)/2,/o} , (4.23) 
{J a } = {-(iV /2-l)/2,...,(iVo/2-l)/2} , 

where |7 | > (N — l)/2. This is a two-parameter family of excited states, called 
(somewhat misleadingly) "particle-hole excitation" by Coll. To understand the ex- 
cited states, we shall in the following consider systems of finite size, rather than taking 
the thermodynamic limit directly. It should however be quite evident how spectra like 
those of fig. 5 develop into a true continuum in the thermodynamic limit (compare in 
particular figs. 6 and 7). In fig. 5 we show numerical results for the energy-momentum 
spectrum of the states (4.23) for a chain of 40 sites. The same states with k — > —k 
are obtained using negative 7 . One notices a sharp minimum in the exctation energy 
at k/ir = 1.1 = 4:kp (this is why I call these excitations "4&p" singlts). In the ther- 
modynamic limit, the gap at Akp vanishes. These excitations are at the origin of the 
power-law behaviour in the density-density correlation function (3.17) around 4&p. 
Moreover, in the bosonization formalism, the Akp density correlations are entirely 
determined by the charge (cf> p ) modes. Consequently, we identify the charge velocity 
Up as the slope of the excitation spectrum of fig. 5 at k = 0. Finally, we notice that 
the overall structure of the spectrum doesn't change much between weak (U/t = 1) 
and strong (U/t = 16) correlation cases. 

"2kp" triplet and singlet states. Excitations of the J's with all A's and Fs real 
are only possible if M < N/2. The simplest excitations of this type are obtained 
considering M = N/2 — 1 which has total spin S = 1 (triplet). Now the restriction 
(4.19) allows N/2 + 1 different J's, i.e. we have two free parameters (the "holes" in 
the J-sequence), leading to sequences of quantum numbers of the form 

{I 3 } = {_iV /2 + l,..., iVo/2} , 
,h = -N /4 + S aul , (4.24) 
J a = J a -i + 1 + 6 a)ai + 6 a)a2 (a = 2,...,M) , 

where 1 < a\ < a 2 < M , and 6 a ^ is the usual Kronecker symbol. Numerical results 
for this type of excitations are shown in fig. 6. Corresponding states with negative k 
are obtained shifting the {7,} in (4.24) by one unit to the left. There now is a sharp 
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Figure 5: %kp" singlet excitation spectrum for a Hubbard chain of J f sites with 22 
electrons. The lowest "arc" from k/ir = 0.05 to k/ir = 1.1 is obtained by varying i 
at fixed I = (N + l)/2 (cf. eq. (4.23)), the higher arches correspond to increasing 
I up to (L - I) /2. 
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minimum at k/ir = 0.55 = 2kp, the gap again vanishing in the thermodynamic limit. 
In the long-wavelength limit, these states are the only spin-carrying excitations at 
constant particle number, so that the slope of the spectrum at k = is equal to 
the spin velocity of the bosonized model. The low-energy excitations around 2kp 
are responsible for the spin contribution to the 2k p spin-spin correlations. As in 
fig. 5, the structure of the excitations doesn't change much between weak and strong 
correlations, however, the energy scale does. In fact, the lowering of the energy scale 
in going to strong correlations corresponds to the lowering of the exchange energy 
(ss 4:t 2 /U) in the strong correlation limit. 

The results of fig. 6 show apparent gaps at 2k p and Akp. These are clearly finite 
size effects, as can be seen comparing with fig. 7 which contains results for the triplet 
excitations for the same particle density and interaction strength as in fig. 6, but for 
a chain four times longer. The gaps are now much smaller and obey approximately 
the 1/L scaling expected. One also easily recognizes how the spectra develop into a 
continuum in the thermodynamic limit. 

Together with the triplet excitations (4.24) there are also singlet states (i\7 = 
iV /2), which are obtained by having one pair of complex conjugate A's among the 
solutions to the original equations (4.14, 4.15). The energies of these states are shown 
by triangles in fig. 5. It is remarkable that these states are nearly degenerate with the 
triplet states and in fact become exactly degenerate in the thermodynamic limit. 

The existence of singlets and triplets with the same energy shows that these 
states are in fact the combination of two noninteracting spin- 1/2 objects, commonly 
called spinons. Of course, because of total spin conservation, these objects can only 
be excited in pairs as long as one keeps the total number of particles fixed. 

Added particle. Adding one particle to the 4z/ + 2 ground state and leaving M 
unchanged, the 7's and the J's are HOI. There are now M + I allowed values for the 
M distinct J's, and the low-energy states then are parametrized by 

{Ij} = {-(iV -l)/2,...,(iVo-I)/2,/o} , 
■h = -M/2 + 8 aul , (4.25) 
J a = J a -i + I + 6 a)ai (a = 2,...,M) , 

where |7 | > (N — l)/2, and I < a\ < M. Corresponding spectra for different 
bandfillings are shown in fig. 8 for I > 0. The symmetric spectra with negative k are 
again obtained using I < 0. The state of minimal excitation energy has momentum 
kp, as in the noninteracting case. The shallow arches in fig. 8 correspond to varying 
«i, and are in fact of the same shape as the lowest branch (from k = 0.057T to 
k = 0.557r) in fig. 6, i.e. they correspond to single-spinon excitations. Close to kp 
the energy of theses states varies as u a (k — kp). On the other hand, varying 7 , one 
goes from one arc to the next, and the corresponding excitation energy (the upper 
limit of the quasi-continuum) varies as u p (k — kp). One also sees that going from one 
arc to the next in fig. 8 (i.e. increasing 7 ) the shape of an individual arc is basically 



24 



L=40, N=22, U/t=l 
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Figure 6: "2kp" spin singlet (A) and triplet (o) excitation spectrum for Hubbard 
chains of J f sites with 22 electrons for different interaction strengths. 
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L= 160,N = 82,U/t= 16 




Figure 7: "2kp" triplet excitation spectrum for a Hubbard chain of 160 sites with 82 
electrons. 

unchanged. Varying I corresponds to a variation of the momentum of the added 
particle, and the figure thus shows that the total energy of a state is just the sum of 
the spinon energy and the "charge" energy associated with the added particle. One 
thus sees that charge and spin degrees of freedom do not interact. This is certainly in 
agreement with the predictions of the bosonization formalism, however, the fact that 
spin and charge separate even in highly excited states is special to the Hubbard model 
(a priori, the bosonized theory can only be expected to be an effective low-energy 
theory for the Hubbard and other lattice models). 

Another notable fact in fig. 8 is the number of available states as I is varied: 
for iVo = 14, 22, and 30 there are respectively f 3, 9, and 5 spinon arches. This means 
that without exciting the spins there are 26, 18, and f states for the extra particle 
available (counting states at negative i.e. the I branch stops at k = ir — kp, rather 
than at k = tt as in a noninteracting system. The explanation of this fact is rather 
straightforward for large U when double occupancy of sites is forbidden: in a system 
with L sites and N electrons, there are only L — N sites available at low energies. 
These states then form the "band" in fig. 8. There are of course states involving 
doubly occupied sites, however these are separated from the continuum of fig. 8 by a 
gap (these states are solutions of (4.14, 4.15) with complex k } s [Woynarovich, 1982]). 
This separation of states occurs for any, even very small U . 
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Figure 8: Excitation spectra for one particle added into a Hubbard chain of 40 sites 
with 14, 22, and 30 electrons. The shallow arches correspond to varying a\ (cf. 
eq. (4.25)) at constant I , and I increases from one arc to the next. Zero energy 
corresponds to the N + 1 particle ground state. 
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Added hole. Finally, let us consider states with one hole in the 4z/ + 2 groundstate 
(N = N - 1, M = (N - l)/2). Then both the J's and the J's are integers. The 
energy is minimized chosing consecutive J's between — (N — l)/2 and (A" — l)/2, but 
there are Af + f possibilities for the Af J's. States corresponding to the sequence 

ih) = {-(iV-l)/2,...,(iV-l)/2} , 
Ji = -M/2 + 6 aul , 

J„ = J a _i + l + 5 ai ,„ (a = 2,...,M) (4.26) 
are shown as the lowest arc between & = — 0.257T and & = 0.257T. This is a one-spinon 
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Figure 9: Excitation spectra for an added hole in a Hubbard chain of J f sites with 22 
electrons. The shallow arches correspond to varying a\ (cf. eq. (4-28) at constant ji, 
and ji increases from one arc to the next. 

branch, this state having necessarily S = 1/2, with velocity u a . One can of course 
also create a hole in the sequence of J's. The energy spectrum for the quantum 
numbers 

h = -(N + l)/2 + S jul 

I, = 1^ + 1 + 6^ (j=2,...,N) (4.27) 

■h = -M/2 + 6 aul , 

J a = Ja-1 + l + ^.a (a = 2,...,M) (4.28) 

(where the free parameters obey 1 < ji < N, 1 < a\ < Af) is also shown in fig. 9. 
Similarly to fig. 8, one notices that varying the "charge" quantum number ji one 
creates a branch with velocity u p . 
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4.2.2 Limiting cases 

Strong correlation. For U — > oo the ratio sinkj/U in (4.14,4.15) clearly vanishes, 
howver there is no such restriction on the A's. Introducing the scaled variables A a = 
2X a /U } the "spin equation" (4.15) becomes 

2A„ + i\ N £- A„ - A + i . 

^ " II a a ■ , (4-29) 



2A a -%) ^ = A 1 A a -A i . 

which are the well-known Bethe ansatz equations for the spin-1/2 Heisenberg chain 
[Faddeev and Takhtajan, 1981], i.e. the spin wavefunction of N particles is just that 
of an iV-site (not Z-site!) Heisenberg chain, even when there is less than one particle 
per site. The "charge equation" (4.14) becomes 

M 

e ikjL = & iP H ^ p H = ^[2arctan(2A„) - tt] , (4.30) 

a = l 

where Pjj is the momentum of the spin excitations of the iV-site Heisenberg chain. 
Eq. (4.30) shows that in the strong correlation limit the allowed &'s are quantized in 
units of 2tt / as for spinless noninteracting fermions (each k is occupied oncel). This 
is of course a consequence of strong correlations: double occupancy is forbidden, i.e. 
the strong local repulsion plays the role of a Pauli principle acting even for electrons 
of different spin (this is only true in one dimension). The boundary conditions deter- 
minig the allowed k } s are determined by the spin momentum Pjj. In particular, for 
4z/ + 2 particles, the ground state momentum of the spins is an odd multiple of 7r, so 
that the allowed k } s are HOI multiples of 2tt/L } and consequently the ground state 
is uniquely determined. 

In this limit the total wavefunction separates into a product of a iV-particle 
Slater determinant for spinless fermions, describing the charge motion, with the spin 
wavefunction of the Heisenberg chain. This separation has been used to calculate 
(numerically) a number of correlation function of the one-dimensional Hubbard model 
in the strong correlation limit [Ogata and Shiba, 1990] 

Weak correlation. For U — > the arctan functions in (4.17, 4.18) become sign 
functions, e.g. 

2 arctan[4(A — sin k)/U] — > 7rsign(A — sin k) . (4-31) 
Eq.(4.18) then becomes 

N M 

sign(A„ - sin kj) = 2J a + ^ sign(A„ - \p) 

3=1 13=1 

= 2J a + 2a - M - 1 , (4.32) 

where in the second line we have assumed (without restriction of generality) that 
{J a } and therefore {A„} are increasing sequences. Because of the restriction (4.19), 
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the minimum value of the r.h.s. of (4.32) is —N + 2, i.e. the smallest sin kj is smaller 
than the smallest X a . From (4.32) one immediately finds 

N 

^[sign(A„ + i - sin kj) - sign(A„ - sin kj)] = 2(J a+1 - J a ) + 2 , (4.33) 

i.e. if J„+i — J a = n there are exactly n + 1 solutions kj satisfying X a < sin kj < X a+ i. 
Proceeding in the same way with (4.17), one obtains 

M 

L(k j+1 -kj) = 27r(7j + i-7j) + 7r ^ [sign(A„ -sin k j+1 ) -sign(A„ -sin kj)] ,(4.34) 

a = l 

i.e. if there is one X a satisfying sin kj < X a < sin k 1+ i one has 

k i+i ~ k i = 2vr(7 J+1 - Ij - 1)/L , 

whereas if there is no such X a 

k j+i ~ k j = 2tt(7 j+1 - Ij)/L . 

From these considerations, the structure of the solution {%}, {A„} correspond- 
ing to the different states above can be obtained in the limit U — > 0. For the ground 
state, one has a sequence {A„}, and associated with each X a there is a pair (k 2a -i } k 2a ) } 
so that 

(sm&2 Q ._i,sin&2 Q .) — > X a =F e , (4.35) 

and e — >■ as U ^0 (see fig. 10). 

The U — > limit of excited states can be discussed similarly. For the Akp 
singlets above there are two adjacent k } s missing in the ground state configuration, 
one of which is put into the first free state above kp, whereas the other goes into a 
state of higher momentum, depending on 7 in (4.23). Obviously, for £7 = these 
states do not contribute to the density correlation function, however, for any nonzero 
U these two particle/two hole states do couple to the particle-hole excitations, and 
therefore these states give rise to the Akp power law in the density correlation function 
(3-17). 

To create the 2k p triplets (4.25) one takes two particles at distinct k } s out of 
the ground state distribution, and puts both into the first available state just above 
kp. The two unpaired electrons then form a spin triplet, or a singlet if one creates a 
2kp singlet. The 7/ — > limit of other excited states can be found in the same way. It 
is important to note here that the simple excited state of the Bethe ansatz discussed 
above correspond in fact to relatively complicated two particle/two hole states of the 
noninteracting system. 
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Figure 10: Dependence of the k's in the ground-state on U. There is one X between 
each pair that converges to a doubly occupied state as U — > 0. 

4.3 Low energy properties of the Hubbard model 

4.3.1 Luttinger liquid parameters 

In a weakly interacting system the coefficients K p and u v can be determined pertur- 
batively. For example, for the Hubbard model one finds 

K p = 1 - U/(7rv F ) + ... , (4.36) 

where vp = 2tsin(7rn/2) is the Fermi velocity for n particles per site. For larger U 
higher operators appear in the continuum Hamiltonian (3.3), e.g. higher derivatives 
of the fields or cosines of multiples of v^V- These operators are irrelevant, i.e. they 
renormalize to zero and do not qualitatively change the long-distance properties, 
but they do lead to nontrivial corrections to the coefficients u vl K p . In principle 
these corrections can be treated order by order in perturbation theory. However, 
this approach is obviously unpractical for large U, and moreover it is likely that 
perturbation theory is not convergent. To obtain the physical properties for arbitrary 
U a different approach is therefore necessary. 

I note two points: (i) in the small-// perturbative regime, interactions renor- 
malize to the weak-coupling fixed point = 0,K* = 1; (ii) the exact solution [Lieb 
and Wu, 1968] does not show any singular behaviour at nonzero U, i.e. large U and 
small U are the same phase of the model, so that the long-range behaviour even of the 
large U case is determined by the fixed point = 0. Thus, the low energy properties 
of the model are still determined by the three parameters u Pj(7 and K p . 

The velocities be obtained from the long wavelength limit of the "4&p" 

and u 2kp v excitations discussed above. In the thermodynamic limit the corresponding 
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Figure 11: The charge and spin velocities u p (full line) and u a (dashed line) for the 
Hubbard model, as a function of the band filling for different values of U/t: for u a 
U/t = 1,2,4,8,16 from top to bottom, for u p U/t = 16,8,4,2,1 from top to bottom 
in the left part of the figure. 
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Figure 12: The Wilson ratio Rw for the one-dimensional Hubbard model, as a func- 
tion of the band filling for different values of U/t (U/t = 16,8,4,2,1 for the top to 
bottom curves). 
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excitation energies are easily found from the numerical solution of a linear integral 
equation [Coll, 1974]. Results are shown in fig. 11 for various values of U/t. Note 
that for U = one has u p = u a = 2tsin(7m/2), whereas for U — > oo u p = 2tsin(7m), 
u<j = (2Trt 2 /U)(l — sin(27rn)/(27rn)). In the noninteracting case u a oc n for small 
n, but for any positive U u a oc n 2 . The Wilson ratio, eq.(3.12) obtained from the 
velocities is shown in fig. 1 1 . For U = one has Rw = 1, whereas for U — > oo Rw = 2 
for ti / 1. 

To obtain the parameter K p from the exact solution note that the gradient of 
the phase field cf> p is proportional to the particle density, and in particular a constant 
slope of <f) p represents a change of total particle number. Consequently, the coefficient 
u p / K p in eq. (3.4) is proportional to the variation of the ground state energy E with 
particle number: 

1 d 2 E (n) = ttu p 

L dn 2 2K P ' { ' ) 

Note that this quantity is the inverse of the compressibility. Equation (4.37) now 
allows the direct determination of K p : E (n) can be obtained solving (numerically) 
Lieb and Wu's [Lieb and Wu, 1968] integral equation, and u p is already known. 
The results for K p as a function of particle density are shown in fig. 13 for different 
values of U/t. For small U one finds in all greement with the perturbative 

expression, eq. (4.36), whereas for large U K p — > 1/2. The limiting behaviour 
for large U can be understood noting that for U = oo the charge dynamics of the 
system can be described by noninteracting spinless fermions (the hard-core constraint 
then is satisfied by the Pauli principle) with kp replaced by 2kp- Consequently one 
finds a contribution proportional to cos(4:kFx)x~ 2 in the density-density correlation 
function, which from eq. (3.17) implies K p = 1/2. One then finds an asymptotic 
decay like cos(2kpx)x~' i l 2 ln 1 ^ 2 (x) for the spin-spin correlations, eq.(3.18), and an 
exponent a = 1/8 in the momentum distribution function. The result a = 1/8 has 
also been found by Anderson and Ren (preprint), and by Parola and Sorella [Parola 
and Sorella, 1990]. Ogata and Shiba's numerical results [Ogata and Shiba, 1990] are 
quite close to these exact values. 

As is apparent from fig. 13, the strong-coupling value K p = 1/2 is also reached in 
the limits n — > 0, 1 for any positive U . For n — > this behaviour is easily understood: 
the effective interaction parameter is U/vp } but vp goes to zero in the low-density 
limit (corresponding to the diverging density of states). The limit n — > 1 is more 
subtle: in this case nearly every site is singly occupied, with a very low density of 
holes. The only important interaction then is the short range repulsion between holes, 
which can be approximated by treating the holes as a gas of spinless noninteracting 
fermions. Using (4.37), one then again finds K p = 1/2. 

We note that in the whole parameter region, as long as the interaction is repul- 
sive one always has K p < 1, which means that magnetic fluctuations are enhanced 
over the noninteracting case. On the other hand, superconducting pairing is always 
suppressed. 
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Figure 13: The correlation exponent K p as a function of the bandfilling n for different 
values ofU (U/t = 1, 2, 4, 8, 16 for the top to bottom curves). Note the rapid variation 
near n = 1 for small U . 
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It should be emphasized here that the results of fig. 13 are valid for n — > 1, but 
not for n = 1. In fact, in that latter case, there is a gap in the charge excitation 
spectrum, as expected from the umklapp term (3.7), and the correlations of cf> p become 
long ranged. Close to half-filling, the asymptotic behaviour of the charge part of 
correlation functions like (3.17) is essentially determined by the motion of the holes. 
Writing the density of holes as p = 1/(1 — n) one then expects a crossover of the form 
[Schulz, 1980] 

(n(x)n(O)) « cos(2A; F x)[l + {pxf^^x' 1 l n - 3/2 (x) (4.38) 

for the 2kp part of the density correlation function, and similarly for other correlation 
functions. Clearly, only for x ^> 1/p are the asymptotic power laws valid, whereas 
at intermediate distances 1 < a; < 1//) one has effectively K p = 0. Clearly, the form 
(4.38) provides a smooth crossover as n — > 1. 

Results equivalent to the present ones can be obtained using the conformal 
invariance of the Hubbard model [Frahm and Korepin, 1990, Kawakami and Yang, 
1990a]. These results have subsequently be generalized to the case with an applied 
magnetic field [Frahm and Korepin, 1991]. 

4.3.2 Transport properties 

The exact solution of Lieb and Wu can also be combined with the long-wavelength 
effective Hamiltonian (3.3) to obtain some information on the frequency-dependent 
conductivity (t(lo). On the one hand, from eq. (3.13) there is a delta function peak 
at zero frequency of weight 2K p u p . On the other hand, the total oscillator strength 
is proportional to the kinetic energy [Baeriswyl et al., 1986]: 

/oo 
a(u)du = -Tv(H kin )/L . (4.39) 
-oo 

Thus, both the weight of the dc peak and the relative weight of the dc peak in 
the total conductivity can be obtained and are plotted in fig. 14. As expected, far 
from half-filling, all the weight in cr tot is in the dc peak. For exactly half-filling the 
dc conductivity vanishes, due to the existence of a gap A c for charge excitations 
created by umklapp scattering, and all the weight is at u > A c . Fig. 2 then shows 
that as n — > 1 umklapp scattering progressively transfers weight from zero to high 
frequency. The crossover is very sharp for small or large U, but rather smooth in 
intermediate cases (U/t 16). This nonmonotonic behaviour as a function of U can 
be understood noting that initially with increasing U umklapp scattering plays an 
increasingly important role. Beyond U /t 16, however, the spinless-fermion picture 
becomes more and more appropriate, and at U = oo one again has all the weight in 
the dc peak. The linear vanishing of <t as n — > 1 implies a linear variation of the 
ratio n/m* with "doping". 

An interesting question is the sign of the charge carriers, especially close to 
the metal-insulator transition. The standard way to determine this, the sign of 
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Figure 14: Top: The weight of the dc peak in (t(lo) as a function of bandfilling for 
different values ofU/t (U/t = 1,2,4,8,16 for the top to bottom curves). 
Bottom: Variation of the relative weight of the dc peak in the total conductivity 
oscillator strength as a function of the bandfilling n for different values of U : U/t = 1 
(full line), 4 (dashed), 16 (dash-dotted), 64 (dotted), and 256 (dash-double-dotted). 
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the Hall constant, is useless in a one-dimensional system. As an alternative, the 
thermopower can be used which is negative (positive) for electron (hole) conduction. 
In general, calculation of the thermopower is a nontrivial task, as the curvature of 
the bands plays an important role, and the approximate form of the Hamiltonian 
(3.3) is therefore insufficient. Moreover, both charge ans spin entropies can play a 
role. However, close to the metal-insulator transition u p <C u a , and therefore the 
entropy of the charge degrees of freedom is much bigger than the spin entropy. In the 
presence of umklapp scattering, which becomes important close to half-filling, the 
charge part of the Hamiltonian can be transformed into a model of massive fermions, 
with energy-momentum relation Sk = i(v 2 k 2 + A 2 ) 1 / 2 [Emery et al., 1976]. A is the 
charge excitation gap created by umklapp scattering. In general, the quasi-particles 
interact, however close to half-filling this interaction can be eliminated [Schulz, 1980]. 
At half-filling all negative energy states are filled, all positive energy states are empty. 
Doping with a concentration n* of holes, some of the negative energy states become 
empty and only states with \k\ > k* F oc n* are filled. Because of the vanishing 
interaction, a standard formula for the thermopower can be used [Chaikin et al., 
1976] and gives 



ri w 2 klT A 2 , , 

S = A 40) 

6|e| v 2 (k* F ) 2 (v 2 (k F ) 2 + A 2 ) 1 / 2 ' K ' ' 

i.e. approaching the metal-insulator transition from n < 1, the thermopower is hole- 
like, whereas obviously far from the transition (ra <C 1) it is electron-like. The exactly 
opposite behaviour can be found for n > 1. 



4.3.3 Spin— charge separation 

The Hubbard model also provides a rather straightforward interpretation of the spin- 
charge separation discussed above. Consider a piece of a Hubbard chain with a half- 
filled band. Then for strong U there will be no doubly-occupied sites, and because 
of the strong short-range antiferromagnetic order the typical local configuration will 
be 

■■■nnnnnn--- 

Introducing a hole will lead to 

■■■nmonnn--- 

and after moving the hole one has (note that the kinetic term in the Hamiltonian 
does not flip spins) 

-uoumuu- 

Now the hole is surrounded by one up and one down spin, whereas somewhere else 
there are two adjacent up spins. Finally, a few exchange spin processes lead to 

-uoutmm- 
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Note that the original configuration, a hole surrounded by two up spins has split into 
a hole surrounded by antiferromagnetically aligned spins ("holon") and a domain- 
wall like configuration, two adjacent up spins, which contain an excess spin 1/2 with 
respect to the initial antiferromagnet ("spinon"). The exact solution by Lieb and Wu 
contains two types of quantum numbers which can be associated with the dynamics 
of the spinons and holons, respectively. We thus notice that spinons and holons 
[Kivelson et al., 1987, Zou and Anderson, 1988] have a well-defined meaning in the 
present one-dimensional case. 

The above pictures suggest that, as far as charge motion is concerned, the Hub- 
bard model away from half-filling can be considered as a one-dimensional harmonic 
solid, the motion of the holes providing for an effective elastic coupling between ad- 
jacent electrons. This picture has been shown to lead to the correct long-distance 
correlation functions for spinless fermions [Haldane, 1981b, Emery, 1987]. For the 
case with spin, this suggests that one can consider the system as a harmonic solid 
with a spin at each site of the elastic lattice (lattice site = electron in this picture). 

Let us now show that this gives indeed the correct spin correlation functions. 
In a continuum approximation, the spin density then becomes 

(t(x) = J2 S mH X ~ X m) , (4.41) 

m 

where the sum is over all electrons. After a Fourier transformation of the delta 
function the spin-spin correlation function becomes 

{<r(x) • <t(0)} = -i- / dqdq' £ e~^(S m • S m , e ^W^)} . ( 4 . 42 ) 

The exchange energy between adjacent spins is always antiferromagnetic, whether 
there is a hole between them or not, and consequently the low-energy spin dynamics 
always is that of an antiferromagnetic chain of localized spins. Under the additional 
assumption that the spin-spin correlations on an elastic lattice depend mainly on the 
average exchange constant and not so much on the fluctuations induced by motion of 
the electrons, the average in (4.42) factorizes into separate spin and charge factors. 
Following the hypothesis about harmonic motion of the electrons, we write x m = 
R m + u m , where R m = m/\l — n | is the average position of the rath electron and u m the 
displacement with respect to this position. Note that the "harmonic solid hypothesis" 
implies that u m+ i — u m is small, but not necessarily u m and u m+ i separately. In the 
averages over atomic positions now all terms with q ^ q' vanish, and one has 

(a(x) ■ <t(0)} « J dqJ2 e~ tqx (S m ■ 5 ml )e ,?(,i '"" ,i '»' ) (e ,?(% ""-')) . (4.43) 

m,m r 

The average over u m in (4.43) has a power law behaviour: 
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with ct(q) oc q 2 } i.e. it has a smooth ^-dependence. On the other hand, the long- 
distance behaviour of the spin-spin correlations of an antiferromagnetic spin chain is 
[Luther and Peschel, 1975, Affleck et al., 1989] 

(S m ■ S m ,) « (-1)™' \ m - m'l 1 ffl 1 / 2 \ m - m '\ . 

Therefore in (4.43) the ^-integration is dominated by terms with q k> iril — n) = 2kp. 
Replacing the weakly ^-dependent exponent ct(q) by a{2kp) one obtains 

(a(x) ■ <r(0)} « J dqe~ iqx ]T (S m ■ S ro /)e , ' 9(iim - R '» ,) |m - m'\ a{2kF) (4.44) 

m,m' 

= cos(2k F x)x- 1 - a{2kF hn 1/2 (x) . (4.45) 

With the identification a(2kp) = K p , this is precisely the result (3.18). We thus have 
shown that the spin-spin correlations of a correlated electron system can in fact be 
understood as those of an elastic lattice of spins. In that picture, the motion of the 
holes then only provides the effective elasticity for the lattice. 

4.3.4 The metal— insulator transition 

It seems worthwhile here to compare the metal-insulator transition in the one-dimen- 
sional Hubbard model with other scenarios for strongly correlated fermion systems 
in higher dimension (see the review by Vollhardt [Vollhardt, 1984]). In the "nearly 
localized" picture, effective mass effects predominate and enhance both the specific 
heat and the spin susceptibility. Consequently, the Wilson ratio (1/(1 + -^0") m Fermi 
liquid language) remains nonzero as the metal-insulator is approached. On the other 
hand, in the "nearly ferromagnetic" (or paramagnon) picture, only the spin suscep- 
tibility is enhanced significantly, and therefore Rw can be much larger than unity. 
The behaviour found here in the one-dimensional case is quite different from both 
these scenarios: generally Rw < 2, and approaching the metal-insulator transition 
Rw 0. This occurs because generally an enhancement of the mass of the charge 
carriers (i.e. a decrease of u p ) has no influence on the spin degrees of freedom (see 
fig.l). This is rather straightforwardly understood in terms of spin-charge decou- 
pling, as explained in the previous section: charge and spin excitations move nearly 
independently of each other, and in particular the spin dynamics is determined by 
antiferromagnetic nearest-neighbor exchange. In particular the spin susceptibility 
remains finite even when the mass of the charge carrier approaches infinity. 

Let us discuss the metal-insulator transition in more detail. The fact that u p 
and cr vanish linearly as n — > 1 seems to be consistent with a divergent effective 
mass at constant carrier density because u p l/m*,cr ~ n/m*. A constant carrier 
density is also consistent with the fact that kp = 7m/2 is independent of U . It is 
not consistent with the hole-like sign of the thermopower as n — > 1 from below, nor 
with the electron-like sign as n — > 1 from above: if the carriers are holes, the carrier 
density is the density of holes: n* = 1 — n. Treating the holes as spinless fermions, 
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as already mentioned before, one expects <t — > because n* — > 0, and 7 — > 00 
because the density of states of a one-dimensional band diverges at the band edges. 
This agrees with what was found explicitly in section 4.3.1. What is not so easily 
understood in this picture is the fact that kp (i.e. the location of the singularity of 
rik) is given by its free-electron value 7m/2, rather then being proportional to n* . 
One should however notice that is given by the single-particle Green's function, 
which contains both charge and spin degrees of freedom. The location of kp then 
may possibly be explained by phase shifts due to holon-spinon interaction. This is 
in fact suggested by the structure of the wavefunction of the exact solution [Ogata 
and Shiba, 1990]. 

The magnetic properties do not agree with what one expects from an effective 
mass diverging as n — > 1: u a and therefor x remain finite. Moreover, the NMR 
relaxation rate would have the behaviour l/7i = aT + f3\/T, where the first (Kor- 
ringa) term comes from fluctuations with q 0, whereas the second term comes from 
antiferromagnetic fluctuations with q 2kp. None of these properties is strongly 
influenced by the diverging effective mass observed e.g. in the specific heat. This 
fact is of course a manifestation of the separation between spin and charge degrees 
of freedom. 

4.3.5 Other models 

For more complicated models, e.g. the "extended Hubbard model" 

H = -t ^2(a] s a i+ltS + 4 +1)S a«) + U ^ n A n ti + V ^ n t n %+1 , (4.46) 

i,s i i 

exact eigenvalues can not be obtained in the thermodynamic limit. The parameters 
in eq. (4.37) can however be calculated reliably for finite systems, and this gives 
rather good results, as shown previously [Schulz, 1990b]. 

Exact exponents can be obtained for the model (4.46) in the limit U — > 00: 
then one has effectively spinless fermions (with kp — > 2kp) with nearest neighbour 
interaction, a model which can be exactly solved using the Jordan- Wigner trans- 
formation into the XXZ spin chain. In particular, the 4^-component of (3.17) is 
related to the correlation function of S z . From the known results [Luther and Peschel, 
1975] one obtains, for a quarter-filled band (n = 1/2), K p = 1/(2 + (4/7r) sin -1 (v)), 
u p = nty/l - u 2 /cos _1 (u), with v = V/2\t\. N ow Kp < 1/2 is possible. For v > 1 
the system is in a dimerized insulating state. Approaching the insulating state from 
v < 1 both K p and u p remain finite, i.e. <t jumps to zero at v = 1. For n / 1/2 
the parameters u p , K p can be obtained from numerical results [Haldane, 1980]. Quite 
generally, one has K p > 1/8, but K p = 1/2 for n — > 0,1, independent of v. On 
the other hand, u p — > as n — > 1/2 for v > 1, i.e. in that case the weight of the 
dc conductivity goes to zero continuously, the point (v,n) = (1,1/2) is thus highly 
singular. The same type of singularity also occurs at U = 0, n = 1 in the Hubbard 
model. Interestingly enough, one has K p > 1 if V < -y/2\t\, i.e. a finite amount of 
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nearest-neighbor attraction is sufficient to lead to divergent superconducting fluctu- 
ations even lor infinite on-site repulsion. Also note that the singularities in u p and 
K p at v = — 1 (attractive interaction) represent a point ol phase separation. 

For the t — J model, there is one exactly solvable point it = J) where exact 
exponents can be found using the Bethe ansatz [Kawakami and Yang, 1990b]. Away 
from this point, eq. (4.37) has been used to obtain K p from numerical data [Ogata 
et al., 1991]. For large J there is a phase with predominantly superconducting fluc- 
tuations. Exact exponents have also been found for a model of fermions with 1/r 2 
interaction [Kawakami and Yang, 1991]. 

5 Conclusion 

In these notes, we have seen that using vastly different techniques, ranging from 
perturbation expansions via the bosonization method to exact solutions, one obtains 
a rather complete picture of many physical properties of interacting one-dimensional 
fermions. Probably the most important feature arising is the Luttinger liquid like 
behaviour, characterized by non-universal power laws, together with the separation 
of the charge and spin dynamics. One should also notice that there are no qualitative 
differences between weak and strong correlation. 

Does this behaviour generalize to higher dimensions? In fact, Anderson has 
suggested that Luttinger liquid behaviour might also occur in two dimensions [An- 
derson, 1990], as well as in coupled chain systems [Anderson, 1991]. At least for 
the coupled-chain case, this is in contradiction with standard scaling [Schulz, 1991] 
and renormalization group arguments [Bourbonnais and Caron, 1991, Fabrizio et al., 
1992]. 

In the weak-coupling limit and close to half-filling, the two-dimensional Hub- 
bard model has some similarities with the one-dimensional case: one has a similar, 
albeit more complicated problem of coupled particle-particle and particle-hole sin- 
gularities [Dzyaloshinskii, 1987, Schulz, 1987]. Upon doping, the antiferromagnetic 
structure becomes an incommensurate spin-density wave [Schulz, 1990c]. However, 
quite unlike the one-dimensional case, upon further doping antiferromagnetism van- 
ishes alltogether. Moreover, strong correlation seems to be quite different from the 
weakly correlated case: for a slightly doped antiferromagnet, both a spiral state 
[Shraiman and Siggia, 1989] and phase separation [Emery et al., 1990] have been 
proposed, but a linearly polarized antiferromagnet seems to be excluded. It thus 
seems that the simple and direct connection between weak and strong correlation of 
the one-dimensional case does not simply carry over to two (or higher) dimensions. 
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